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We consider the slowly rotating relativistic stars with a uniform angular velocity in the scalar- 
tensor gravity, and examine the rotational effect around such compact objects. For this purpose, we 
derive a 2nd order differential equation describing the frame dragging in the scalar-tensor gravity 
and solve it numerically. As a result, we find that the total angular momentum is proportional to the 
angular velocity even in the scalar-tensor gravity. We also show that one can observe the spontaneous 
' scalarization in rotational effects as well as the other stellar properties, if the cosmological value 

of scalar field is zero. On the other hand, if the cosmological value of scalar field is nonzero, 

■ the deviation from the general relativity can be seen in a wide range of the coupling constant. 

■ Additionally, we find that the deviation from the general relativity becomes larger with more massive 
stellar models, which is independent of the cosmological value of scalar field. Thus, via precise 
observations of astronomical phenomena associated with rotating relativistic stars, one may be 

I possible to probe not only the gravitational theory in the strong-field regime, but also the existence 

\ of scalar field. 
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K ■ I. INTRODUCTION 

The general relativity proposed by Einstein is the theory describing the gravity. In order to verify the general 
Q ■ relativity, many experiments have been done up to now and there is nothing indicating the failure of this theory. 
$-H ■ However, the most of these experiments are performed in the weak gravitational field, such as the solar system. 
C/3 , This means that the gravitational theory in the strong gravitational field is still unconstrained by the astronomical 

■ observations. Namely, the gravitational theory might be different from the general relativity in the strong gravitational 
] field. If so, one can probe the gravitational theory by observing the deviation from predictions in the general relativity. 

^-H ■ Actually, the development of technology will enable us to observe the phenomena related to compact objects with high 
I> [ accuracy. Such observations may be used as a test of the gravitational theory in the strong gravitational field. So far, 
. many attempts and possibilities to test the gravitational theory in the strong gravitational field have been suggested 

■ [ir 6] . Hopefully the reliable gravitational theory even in the strong gravitational field will be observationally revealed 
in the near future. 

The scalar-tensor gravity is a natural extension of the general relativity, and one of the simplest alternative gravi- 
tational theory, where the gravity is mediated by long-range scalar fields in addition to the usual tensor field in the 
general relativity [l|, 0| • The scalar field plays an essential role in explaining the accelerated expansion phases of the 
universe, such as inflation scenario. Additionally, the scalar-tensor theory in gravity can be obtained from the low en- 
ergy limit of string and/or other gauge theories. In the scalar-tensor gravity, one adopts the basic assumption that the 
scalar and gravitational fields, ip and 5*^1/, are coupled to matter via an effective metric defined as g^i, = A'^{ip)g^:^^, 
where the function form of the conformal factor is still unclear. The Brans-Dicke theory is the simplest version of the 
scalar-tensor gravity, which assumes that A{ip) — exp(aiy9) The parameter in the conformal factor, a, can be often 
Ci [ connected to the Brans-Dicke parameter, wbd, as = (2wbd + 3)~^, and the expected deviation from the general 
relativity becomes of the order The experiments in the solar system set severe constraints in the parameter, 

i.e., wbd ^ 40000 or a < 10~^ [13] ■ On the other hand, Damour and Esposito-Farese adopted the different function 
form of the conformal factor, i.e., A{(p) = exp{a(p + fiip^ /2). Then, they showed that the massive neutron stars in the 
scalar-tensor gravity can significantly deviate from the expectations in the general relativity [3, [llj . Especially, they 
pointed out that the sudden deviation in the neutron star properties from the predictions in the general relativity 
is possible for the specific values of the coupling constants, even if the value of a is considerably small. This phe- 
nomenon is referred to as "spontaneous scalarization." Subsequently, Harada studied more details with the technique 
of catastrophe theory, and found that the spontaneous scalarization is possible for /3 ^ — 4.35 [T2|. Recently, Freire 
et al. have ruled out the parameter range of /3 < —5 with the observations of the pulsar white dwarf binary [l3|. This 
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constraint is quite strict, but we still have a little chance to see the spontaneous scalarization. So, we focus on the 
range of /3 > — 5 in this article. 

Until now, several possibilities to distinguish the scalar-tensor gravity from the general relativity by using astronom- 
ical observations are suggested, e.g., with the redshift lines of the X- and 7-rays radiated from the surface of neutron 
stars i l4i] and with the spectrum of gravitational waves emitted from neutron stars [l5i . 1I61] . In this article, we examine 
the different approach to probe the gravitational theory. That is, we focus on the rotational effect around neutron 



stars, where we especially consider the slowly rotating neutron stars with a uniform angular velocity 17[. The similar 
analysis in the general relativity has been originally done by Hartle [l9j , and subsequently a lot of studies associated 
with the stellar rotation have been examined. For instance, taking into account the rotational effect, the additional 
oscillation famil y re ferred as r-modes can be excited, and one might be possible to obtain the stellar information from 
their spectrum f20l [2T| . In this sense, this article could become a first step to consider the rotational effect around 
relativistic stars in more complicated system in the scalar-tensor gravity. We should remark that the similar analysis 
in the tensor- vector-scalar (TeVeS) theory has been done ||6|. However, in such an analysis, they did not consider the 
dependence on the scalar coupling, because in TeVeS the dependence on the vector coupling is stronger than that on 
the scalar coupling. In addition, for simplicity, they did not consider the dependence on the cosmological value of 
scalar field. So, in this article, we will see the dependence on the scalar coupling as well as the dependence on the 
cosmological value of scalar field, although the scalar-tensor gravity is one of the specific case in the TeVeS. 

This article is organized as follows. In the next section, we mention the fundamental parts of the scalar-tensor 
gravity and the equilibrium of non-rotating relativistic stars in the scalar-tensor gravity. In Sec. IIIIl we derive the 
differential equation describing the frame dragging in the scalar-tensor gravity. Additionally, we numerically show 
the rotational effect in this section. Then, we make a conclusion in Sec. IIVI In this article, we adopt the unit of 
c — G — 1, where c and G denote the speed of light and the gravitational constant, respectively, and the metric 
signature is {-,+, +, +). 

II. STELLAR MODELS IN SCALAR-TENSOR GRAVITY 
A. Scalar-Tensor Gravity 

The scalar-tensor gravity is a natural extension of the general relativity, where the gravity is mediated not only by 
metric tensor but also by a massless scalar field (p. The total action of such a gravitational theory is given by 01 

S = -7^ / V^{R* - 25r¥',M<^,^') d^x + S„r [*™, A^ip)g,f,,] , (2.1) 



167rG* _ 

where G* is the bare gravitational coupling constant, 5*^1/ is the Einstein metric, i?* is the Ricci scalar constructed 
with represents matter fields collectively, and Sm denotes the matter action, respectively. It is known that 

the formulation of field equations is better in "Einstein frame" described by the Einstein metric. But a test particle 
moves on the geodesic in "Jordan frame" described by the metric 5^^, which is defined as 5^11, = ^'^i'P) g*^Lv |2l- So, 
the Jordan frame is often referred to as physical frame. Hereafter, the quantities with asterisk are related to the 
Einstein frame, while the quantities with tilde denote those in the physical frame. 

By varying the total action S with respect to g*^l,l,^ one can obtain the field equation for the tensor field; 

G^^,y = SttG^T^^^ + 2 (^(p,f^(p^^ - ^g*f,,^g*'^'P,af.,i3^ , (2.2) 

where T^^j^ is the energy-momentum tensor in the Einstein frame, which is associated with that in the physical frame 

as 

rr = ^=T^ = A6((^)f ^^ (2.3) 

Similarly, by varying S with respect to Lp, one obtains the field equation for the scalar field; 

□ *(^ = -47rG*a((/?)r*, (2.4) 
where denotes the d'Alembertian of 9^^^^, and a{p) and are defined as 

, , (ilnA(o?) 

a(^) ^ — (2.5) 



n = T^J- = Tr9*>.u. (2.6) 
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We should remark that the scalar-tensor gravity for a{(p) = reduces to the general relativity. In addition to Eqs. 
(|2.2p and (I2.4p . one gets another equation from the energy- momentum conservation law, i.e., V^T'^ — 0, such as 

V*^T,^, = a((^)r,V,,¥>. (2.7) 

Especially, in this article, we adopt the same form of conformal factor A{(p) as in Ref. @, i.e., 

A{ip) = e^'P^I^, (2.8) 

where /3 is a real constant. Thus, ol{}p) can be written as a((p) = and the scalar-tensor gravity with /3 = 
corresponds to the general relativity. 



B. Nonrotating Relativistic Stellar Models in Scalar- Tensor Gravity 



The equilibrium configurations of non-rotating relativistic stars in the scalar-tensor gravity have been calculated in 



Refs. 



The metric describing the static, spherically symmetric relativistic stars can be written as 



= _e2*('-)di2 + e2A(.)^^2 ^ ^2(^02 ^ gj^a 0^02), (2.9) 

where e"^'^ = 1 — 2\i(r)lr, and /i(r) is corresponding to the mass function. Therefore, the physical metric is 

ds^ ~ Qp^ydx^ dx^ 



A^e^'^dt^ + A'e^'^dr' + A'r'{de' + sin" 0rf0"). 



l2„2A . 2 



(2.10) 

We assume that the stellar matter is a perfect fiuid, i.e., — (^p + Ujj^ + Pgp,u, where p, P, and are the 

total energy density, the pressure, and the four-velocity of the fiuid in physical frame, respectively. 

Using Eqs. (|2.2I) . (|2.4p . and (|2.7p . one can obtain the Tolman-Oppenheimer-Volkoff (TOV) equations in the scalar- 
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(2.12) 
(2.13) 
(2.14) 

(2.15) 



where the prime denotes a derivative with respect to r. In order to close the system of equations, one needs to prepare 
the relationship between p and P, i.e., equation of state (EOS). In this article, we adopt the polytropic EOS as in 
Refs. 0,[il[i5l; 



(2.16) 
(2.17) 



P = Knomii — 
P 

p — nnrib + 



r 



where h is the baryon number density in the physical frame, while m;, and riQ are some constants given by ruh = 
1.66 X 10"^'' g and no = 0.1 fm~'^. In order to fit this EOS to the tabulated data of realistic EOS known as EOS A 
[13 and EOS II [H], we especially adopt T = 2.46 and K = 0.00936 for EOS A and T = 2.34 and K = 0.0195 for 
EOS II. 

The stellar equilibrium can be determined as follows. With the boundary conditions at the stellar center as /i(0) — 0, 
$(0) = $c, p(0) = pc <p(0) = and \I'(0) = 0, the TOV equations with EOS are integrated outward and the stellar 
surface is determined at P = as i? = A{ifs)rs^ where and denote the values of (p and r at P = 0. For r > Ts, 
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the TOV equations with p = P = are integrated outward. And then, the central values of $c and ipc can be 
determined so that the calculated solutions should agree with the asymptotic behaviors as 



2Madm 



o 



(2.18) 
(2.19) 



where rj^^, A/adm, fo, and Q are corresponding to the Minkowskian metric, the ADM mass, the cosmological value 
of scalar field, and the total scalar charge, respectively ^23]. In practice, with the matching conditions at the stellar 
surface between the interior and exterior solutions, one can determine the values of MadMi V'o, and Q as 



Madm = — 1 



(P0 = ips + 



1/2 



exp 



$1 



: arctanh 



arctanh ^ , ^ ; — - 



l/rs 



(2.20) 

(2.21) 
(2.22) 



where the quantities with the subscript s denote their values at r = [o]. At last, choosing the values of two constants 
(3 and (/?o, the stellar models become one parameter family with pc (or Madm)- Especially, in this article, we adopt 
ipo = 0.01 as a example for (pQ ^ 0. 

We show the dependence of the stellar radius R on the coupling constant — /3 with the fixed ADM mass in Fig. 
[T] for EOS A and in Fig. [5] for EOS II. In both figures, the left and right panels are corresponding to the case of 
y^o = and y^o = 0.01, respectively. From these figures, one can observe the spontaneous scalarization in K at around 
/3 ~ —4.4 for — 0, while the stellar radii are changing smoothly for p^ ^ 0. Namely, even for —4.4 < /3 < 0, the 
stellar properties in the scalar-tensor gravity with i^o 7^ are different from those in the general relativity (for the 
case of /3 = 0). 




FIG. 1: For EOS A, the stellar radius R is plotted as a function of — /3 with the fixed ADM mass for <^o = (left panel) and 
for ipo ~ 0.01 (right panel). 



III. SLOWLY ROTATING RELATIVISTIC STARS IN SCALAR-TENSOR GRAVITY 



A. Rotating Dragging 



In order to see the rotational effect around the relativistic stars in the scalar-tensor gravity, we consider a slowly 
rotating stellar model with a uniform angular velocity, fl, where the rotational axis is put on 9 — 0. For this purpose, 
we assume to keep only the linear effects in the angular velocity, which is the same treatment as in the general relativity 
[l9| . With this assumption, the stellar model is still spherically symmetric, because the deformation of stellar shape 
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FIG. 2: Similar to Fig. [T] but for EOS II. The labels denote the adopted stellar mass in the unit of Mq, where the seven solid 
lines from top to bottom correspond to the cases of M/Mq — 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, and 1.9. 



due to the rotation is of the order fj^. Similarly, the rotational effects in the components of metric except for gt^ 
become of the order fl^. So, the metric describing slowly rotating system in the physical frame can be written as 



dS"^ = -A^e^'^dt^ + A^e^^dr^ + A^r^{de^ + sm^ 6 d4)'^) - 2ur^A^ sin^ edtd<j), 



(3.1) 



where the last term appears due to the rotational effect. Hereafter, in order to specify the rotational effects, the 
deviation from the non-rotating stellar model will be expressed by using variables with (5, such as 6gt,p = —ujr^A^ sin^ Q. 
The deviations of the pressure, density, and scalar field due to the rotation are also of the order ©(Jl^), i.e., 5P — 
0(17^), 8p — 0(^2^), and 8^ = 0{Vt^), because those properties should behave in the same way under a reversal in 
the direction of rotation. With such a ordering, one can show that 5g^i, — A^Sg^^^i, and ST^^i, = A^dT^^, within the 
linear order of £7. Additionally, the deviation of fluid velocity in the physical frame can be described as 

5?7^ = (0,0,0,f2?7*) , (3.2) 

where f/* denotes the i-component of four velocity for non-rotating stellar model, i.e., C/* — e^*/A. Consequently, 
the nonzero component of the energy-momentum tensor of the order VL is only (JTt^, which are given as 



5ft^ = r'^A^ ( p + P 



[uj - sin^ d - Pi^r^A^ sin^ Q. 
At last, the rotational effect w can be determined by solving the following equation; 



-2A,T,2 



(3.3) 



(3.4) 



which comes from the Einstein equation (|2.2p . We should remark that another field equation (|2.4p can not tell us any 
information about the rotational effect in the order of (l, because 6ip ~ 0{(l^), 6a ~ 0{d'^), and ST^ ^ 0{Cl^). One 
can show that only the {t, (j)) component in Eq. (j3.4l) becomes a nontrivial equation, which can be written down as 



$' - A' uj' 











j,2 













$" + ($' _ A') 



uj.ee + 'iuj e- 



cos 9 



lQT:G,A^e^^ 



0. 



(3.5) 



In general, uj ~ 0{(l) can be expressed as a;(r, 9) = —uj{r)/ sm9 dgPe{cos9), where Pg is the Legendre polynomial of 
order £. With this expression, the bracket of the fourth term in Eq. p.Sp can be reduced to 



cos 9 

LOfie + oLOfi— — - 



(£ + 2)(£-l)c.55 



(3.6) 



Furthermore, in order to reproduce the general relativity limit of the scalar-tensor gravity, i.e., /? — >■ 0, we especially 
adopt the case for £ = 1 in this article. In practice, one can not adopt an arbitrary value of £ except for ^ = 1 to 
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construct the distribution of w which satisfies the regularity condition at the stehar center [19j . Finally, the equation 
to solve becomes 

J' + (^1 - $' - A'^ J - 2 

Note that Eq. (13.71) for /3 = agrees with the well-known equation describing the frame dragging in the general 
relativity 19]. 

B. Numerical Results 

In order to determine the rotational effect a;(r) with Eq. (13. 7p . one needs to prepare the boundary conditions at 
the stellar center and the spatial infinity, i.e., the regularity condition at r = and the asymptotic flatness far from 
the central object. Since one can expand w as ui(r) = tUco + i^c2t'^ + • • ■ in the vicinity of stellar center, Eq. p.7p is 
integrated outward from r = with the boundary conditions of w(0) = cjco a-nd w'(0) = 0. Then, one can find the 
exact value of lOcq in such a way that the solution uj{r) should satisfy the asymptotic behavior, i.e., 

= (3.8) 

where J is some constant. In practice, we adopt the boundary condition of 1 + Zuj/{ui'r) = 0, which can be derived 
from Eq. p.8p . at the numerical boundary far from the central object. In the general relativity, the distribution of 
ijj{r) outside the star can be analytically written down as a;(r) = 2J/r^, where the constant J corresponds to the 
total angular momentum of the central object, and J is related to the angular velocity for slow rotation as 

J = m. (3.9) 

Here, the constant of proportionality / corresponds to the relativistic generalization of momentum of inertia for slowly 
rotating system (l9j . On the other hand, in the scalar-tensor gravity, we plot the value of J /Vl as a function of the 
angular velocity 17 for the stellar models with EOS A and with A/adm = \AMq in Fig. [3] for ipQ — Q and in Fig. 
|3]for (po 7^ 0. In the both figures, the right endpoints of lines are corresponding to the allowed maximum angular 
velocity for each stellar model. From these figures, one can see that the value of J/O is independent of Q, even in the 
scalar-tensor gravity. That is, Eq. p.9|) can hold even in the scalar-tensor gravity. This is the same results as in the 
TeVeS Q. We remark that in Fig. [3] the lines for /3 > —4.4 are completely equivalent to the line for /3 = 0, because 
the stellar models with iy9o = for /3 > —4.4 are the same as those in the general relativity (for /? = 0). 
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FIG. 3: Dependence of J /Cl on the angular velocity Cl for the stellar models with EOS A and with A/adm = 1.4M0 and tpo — 0, 
where the solid line corresponds to the case in the general relativity, while the broken lines correspond to the results in the 
scalar-tensor gravity with different values of j3. 

The distribution of cj{r) for the stellar model with EOS A and with Madm = 1.4M0 and 17 = 1 kHz is plotted in 
Fig. [5] for ipo — and in Fig. [6] for ipo = 0.01. In the both cases, we can obviously observe the difference of uj inside 
the star in the scalar-tensor gravity from the prediction in the general relativity. In fact, compared with the case in 
the general relativity, one can see the decrease of the central value of lo in the scalar-tensor gravity up to 6.3% in Fig. 
[5] and 10.0% in Fig. [S] Thus, via the observations such as the stellar oscillations and/or the radiated gravitational 
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FIG. 4: Similar to Fig. [3l but for ipo = 0.01. The labels attached on each line denote values of p. 



waves, one could be possible to probe the gravitational theory in the strong-field regime. Meanwhile, although the 
deviation outside the star may be almost negligible, there still exists. Now, we should remark the difi^erence between 
the results in the scalar-tensor gravity and in the TeVeS. Unlike the results in the scalar-tensor gravity, the value of 
uj{r) in the TeVeS becomes larger than that in the general relativity in the vicinity of stellar center, and becomes 
smaller in the outer region (see Fig. 3 in [^). This means that one might be possible to probe the gravitational theory 
via the observations of binary system composed of the relativistic stars with quite high accuracy. 
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FIG. 5: Distribution of cj(r) for = with different values of /3, where the adopted stellar model is A/adm = 1.4M0 and 
j7 = 1 kHz for EOS A. 




FIG. 6: Similar to Fig. [5] but for ifio = 0.01. 



Furthermore, we plot the momentum of inertia / as a function of the coupling constant —f3 with the fixed ADM 
mass in Fig. [7] for EOS A and in Fig. |8] for EOS II, where the left and right panels are results for ipo = and for 
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ipo = 0.01. From the both figures, it is found that the dependence of / on — /3 becomes stronger, as the stellar mass 
increases. One can also observe that, compared with the case of (/Jq =0, the case oi ip^ ^ Q depends strongly on the 
scalar coupling /3. For example, for EOS A, the value of / in the scalar-tensor gravity becomes up to 1.4%, 5.3%, 
11.0%, and 19.7% larger than that in the general relativity for M = l.SM©, 1.4Mq, I.SMq, and 1.6Mq in the left 
panel of Fig. H while up to 4.7%, 9.9%, 16.8%, and 27.4% larger for M = 1.3Mq, 1.4Mq, 1.5Mq, and 1.6Mq in the 
right panel of Fig. [71 We should emphasize that such a dependence of / on the gravitational theory is completely 
different from that in the case of TeVeS, i.e., the value of / in the TeVeS becomes smaller than that in the general 
relativity Additionally, for the case oi (po = (left panels in Figs. [7]and|S]), only if /3 is less than around —4.4, / 
in the scalar-tensor gravity can deviate from / in the general relativity. Namely, the spontaneous scalarization can be 
observed in / as well as the other stellar properties. On the other hand, for the case of (^o (right panels in Figs. 
[7] and [5]), one can observe the deviation in I even for —4.4 < /3 < 0. Now, in order to clarify the deviation in / for 
(fio = 0.01, wc define the relative deviation A as 

A - {1st - /gr)//gr, (3.10) 

where 1st and Iqr denote the momentum of inertia defined as Eq. p.9|) in the scalar-tensor gravity and in the general 
relativity, respectively. The calculted A is plotted as a function of — /? in Fig. ^ where the left and right panels are 
corresponding to the results for EOS A and for EOS II. This figure obviously shows the statement mentioned the 
above, i.e., the deviation from the general relativity becomes larger with more massive stellar model. Meanwhile, 
comparing the stellar models with same mass but different EOS, one can see that the stellar model with softer EOS 
becomes larger deviation from the prediction in the general relativity. At last, one can clearly see the deviation from 
the general relativity even for —4.4 ^ /? < 0, which is around 0.01 — 10 %. This deviation might be small, but we still 
have a chance to probe the gravitational theory in the strong-field regime via the accurate observations around the 
relativistic objects. 




FIG. 7: For EOS A, the momentum of inertia / is plotted as a function of — /3 with the fixed ADM mass for ipo — (left panel) 
and for ifo — 0.01 (right panel). 




FIG. 8: Similar to Fig. [3 but for EOS II. 
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FIG. 9: Relative deviation of / in the scalar-tensor gravity compared with I in the general relativity. The left and right panels 
correspond to the case of EOS A and EOS II, respectively. 



IV. CONCLUSION 

In order to examine the rotational effect around the relativistic stars in the scalar-tensor gravity, we consider 
the slowly rotating relativistic objects with a uniform angular velocity. From the Einstein equations, we derive a 
second order differential equation governing the frame dragging, and solve it numerically with appropriate boundary 
conditions. As a result, we find that, similarly in the general relativity, the total angular momentum is proportional 
to the angular velocity for the slowly rotation. Additionally, we find that the spontaneous scalarization can arise 
in the rotational effects as well as the other stellar properties for (^sq = 0, where lpq denotes the cosmological value 
of scalar field. On the other hand, for (^o 7^ 0, we can observe the deviation from the general relativity in a wide 
range of the coupling constant /3. Comparing the results in this article with those in the case of TeVeS, we find 
that the dependence on the gravitational theory is completely different. Especially, we find the obvious difference in 
the momentum of inertia. That is, compared with the prediction in the general relativity, the momentum of inertia 
in the TeVeS becomes smaller, while that in the scalar-tensor gravity becomes larger. Thus, via the astronomical 
observations around the relativistic stars with high accuracy, one could possible to probe not only the gravitational 
theory in the strong-field regime, but also the existence of scalar field. At last, we adopt simple stellar models in 
this article, but to compare with the actual observational data, we need to consider the more realistic stellar models, 
where we might take into account the magnetic effects [1^ and the effects of crust and/or the more exotic phase 
[i^ll^]. Considering such additional properties, one could be possible to obtain the further constraint in the theory. 
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